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SUMMARY 


A series representation of the oscillatory behavior of incompressible non- 
viscous liquids contained in partially filled elastic tanks is presented- By 
selecting each term of the series on the basis of the hydroelastic vibrations in 
circular cylindrical tanks, each term satisfies the governing liquid equation 
(Laplace's equation) but does not satisfy the liquid-tank interface compati- 
bility. By using a complementary energy principle presented herein, the super- 
position of these terms is made to approximately satisfy the interface compati- 
bility. This analysis is applied to the gravity sloshing and hydroelastic 
vibrations of liquids in hemispherical tanks and in a typical elastic aerospace 
propellant tank. With only a few series terms being retained, the results 
correlated very well with existing analytical results, NASTRAN^-generated 
analytical results, and experimental test results. Hence, although each term 
is based on a cylindrical tank geometry, the superposition can be successfully 
applied to noncylindrical tanks . 


INTRODUCTION 

Liquid-fuel tanks are often large components of aerospace vehicles and 
significantly affect the vibratory behavior of an entire vehicle. In turn, the 
motions of liquid in a fuel tank can be expected to play a significant role in 
tank vibrations. Consequently, vibration analyses of such tanks must include 
interactive coupling between the elastic shell tank structure and the contained 
liquid. 

Early analytical investigations of liquid-tank vibrations were mostly 
limited to simple tank geometries, such as circular cylinders. A thorough sur- 
vey of these early investigations is given in reference 1, chapter 9. More 
recently, computer programs have been developed which treat the liquid- tank 
interaction problem in more general tank geometries found in practice. (See, 
for example, refs. 2 to 7.) The procedure used in these references to develop 
such programs was to modify existing shell programs to include liquid- structure 
interaction. The key element in these interaction analyses is the selection of 
a suitable representation of the liquid oscillatory behavior. 

Since modeling the shell tank alone may require many degrees of freedom, a 
liquid model which introduces a minimum number of additional degrees of freedom 
is desirable. In references 2 to 5, finite-element models of the liquid are 
employed. Such models introduce a large number of degrees of freedom into the 
interaction problem since, in general, a three-dimensional body of liquid must 
be modeled. Furthermore, the distribution of finite elements throughout the 
liquid may pose modeling problems. In reference 6, the liquid is represented 
by a distribution of simple sources on the surface of the liquid. With this 
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approach, only the liquid surface is discretized and, consequently, the number 
of degrees of freedom is drastically reduced from the number introduced by the 
finite-element approach. In reference 7, the liquid model involves a series 
representation of the liquid oscillatory pressure with the liquid considered 
incompressible and gravity surface, effects neglected. In this model, the 
degrees of freedom are the amplitude coefficients of terms of the series. The 
series is chosen so that each term exactly satisfies the governing equation for 
the behavior of the liquid (Laplace's equation). Using a series composed of 
these terms in a complementary energy principle for the liquid produces surface- 
integral governing equations. Thus, as in reference 6, only the liquid surface 
is discretized for numerical solution of these equations. Furthermore, the 
individual terms of the series capture the essence of the liquid behavior; con- 
sequently, the total number of degrees of freedom required for the coupled 
liquid-shell interaction problem is only slightly greater than that required 
for the empty shell- In addition, this approach can be applied successfully to 
a wide variety of shell geometrical shapes. 

Inasmuch as the pressure-series representation provides accurate solutions 
when gravity effects are neglected, it is reasonable to anticipate that the 
approach can be successfully extended to include such effects. The inclusion 
of gravity permits the analysis to be applied to the physically important slosh- 
ing modes of the liquid which occur with little or no interaction with the 
elastic tank. 

The purpose of the present paper is to extend the analysis of reference 7 
to include gravity effects in the liquid. Further, since the analysis and 
application to incompressible liquids were only briefly described in reference 7, 
a more thorough description is presented herein. A series representation of the 
liquid oscillatory pressure is employed in a complementary energy principle 
which is derived herein for incompressible liquids in the presence of gravity. 

The terms of the series are uniquely selected to satisfy the governing liquid 
equation. Then, as in reference 7, the liquid pressure loading on the tank wall 
which depends on the wall deformation may be calculated without further dis- 
cretization of the liquid. The deformation-dependent pressure of the liquid 
on the tank wall then becomes a loading term in the shell equations, and existing 
shell computer programs may be easily modified to accommodate this loading. In 
the present paper, the modification of the shell-of-revolution program of ref- 
erences 8 and 10 is discussed. 

Applications of the modified shell program are presented for the following 
problems: liquid sloshing in a rigid hemisphere, hydroelastic vibrations of an 

incompressible liquid in an elastic propellant tank, and hydroelastic vibrations 
of an incompressible liquid in an elastic hemisphere. Furthenmore , the conver- 
gence of the series expansion is considered in each application, and the results 
are compared with experiment and with NASTRAN analyses to validate the approach. 
Appendixes to the report provide a proof for the complementary energy principle 
used herein and descriptions of the NASTRAN analysis and experiment. 
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B,B* ,B*' 
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9f9 

h 

j 

Jn 
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m 

M 

' ^ 2 ' ^3 


n 


radius of hemisphere 

A matrices whose elements are defined by equations (32) to (35) 

mth amplitude coefficient of series expansion, see equations (8) 
and (9) 

matrices whose elements are defined in equations (51), (41), and (46), 

respectively 

defined by equation (19) 

diagonal matrix whose elements are 

column vector of forces from liquid acting on the wetted shell surface 
at l£ nodal stations 

gravitational vector in x-direction and magnitude, respectively 
height of liquid 

number of nodal stations on wetted shell surface at which liquid 
pressures and forces are calculated 

modified Bessel function of first kind of order n 

= nPT 

Bessel function of first kind of order n 
modified Bessel function of second kind of order n 
length of tank in axial direction 
number of half-waves along meridian 

half the total number of terms in series expansion 
liquid apparent mass matrix 

number of series terms associated with positive, negative, and zero 
values of respectively 

number of circumferential waves 
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So 

Sf 
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T 


c 


u 
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c,g 
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^mn 


component of N in x-direction 

outward unit vector normal to shell surface 

liquid oscillatory compressive pressure 

mth series term in nth circumferential trigonometric harmonic 
series term defined in equation (61) 

matrices whose elements are found from equation (34) 
radial coordinate shown in figure 1 
value of r at shell wall 

meridional coordinate of shell shown in figure 1 
shell surface area 

free and wetted shell surface, respectively, of liquid 
time 

transformation matrix defined by equation (51) 
complementary kinetic energy of liquid 
liquid displacement vector 
component of u in x-direction 
gravitational energy 
liquid volume 

outward normal displacement of shell 
complementary work of applied displacements 
axial coordinate shown in figure 1 
position vector of liquid particle, equation (7) 

Bessel function of second kind of order n 
defined in equation (17) 
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Tm ' Itt 


■m 
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defined by equations (15) and (16) , respectively 
gravitational load matrix defined by equation (53) 
variational operator 

Kronecker delta function defined in equation (35) 

= g/oj^h 

eigenvalue of equation (46) 

circumferential coordinate shown in figure 1 

mth value of selected parameter which determines type of series term 
to be used according to equations (14) to (17) 

mass density of liquid 

defined by equation (20) 

velocity potential for liquid 

liquid load matrix defined by equation (52) 

angle between N and r axis shown in figure 1 

circular frequency 

derivative on shell wall parallel to N 
A dot over a symbol denotes differentiation with respect to time. 

A prime denotes differentiation with respect to x. 

Bars over symbols refer to associated dimensionless quantities given by 
equation (27) . 
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GENERAL THEORY FOR OSCILLATING LIQUID 
Variational Principle 

In appendix A, the following complementary energy principle is presented 
for a nonviscous incompressible liquid undergoing irrotational oscillations of 
the complex form e^^^ about an initial compressed state: 
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where the variations indicated by the operator 6 are on all admissible pres- 
sure states, 0) is the oscillatory circular frequency, j is and the 

variations of the complementary kinetic, gravitational, and external work 
energies are given by 
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dS, 
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In equations (2) to (4) , p is the oscillatory compressive pressure perturba- 
tion above the initial pressurized state P^; is the uniform mass density 

of the liquid; and are the liquid surface area and liquid volume, 

respectively; N is a vector normal to Sq and positive outward; is the 

component of N in the positive x-direction (see fig. 1); g is the gravity 
vector, which acts in the x-direction; refers to the wetted portion of S^; 

w is the shell displacement perturbation on parallel to N; and (j) is a 

velocity potential for the irrotational nonviscous liquid defined such that 


u = grad cj) 


(5) 


where the dot represents differentiation with respect to time and u is the 
liquid displacement vector. 

The admissible pressure states for equation (1) must satisfy the free sur- 
face condition, namely = 0 on S^, the free surface portion of S^; and the 

dynamic equilibrium of the liquid (or equivalently the linearized Bernoulli 
equation as established in appendix A) , namely 


p = p (-4) + 5 
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( 6 ) 


Series Expansion 

The velocity potential is approximated by a truncated series: 
2M 

^ t>m 

m=l 


4>(Y) = jtO 


(7) 



where Y is the position vector of a liquid particle. In equation (7) and 
henceforth, the complex multiplier e^^^ is understood. Substituting equa- 
tions (5) and (7) into equation (6) yields the series expansion for p, 

2M 

p(Y) = Po ^ 
m=l 
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qm(Y) + g 




( 8 ) 


where the prime denotes differentiation with respect to x. Equations (7) 
and (8) may now be substituted into equations (2) to (4) and the results substi- 
tuted into equation (1) to yield 
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( 9 ) 


where the effects of ullage volume have been neglected and the integral over 
vanishes, since 6p vanishes on Sf. 

Since the numerical evaluation of the volume integral occurring in equa- 
tion (9) can be a tedious task, it is desirable to cause this integral to 
vanish. This can be done by judiciously choosing q^^. 


Selection of Series Terms 


From the derivation of the complementary energy principle presented in 
appendix A, the volume integrand is equivalent to a statement of compatibility 
or incompressibility, namely 

V^cj) = 0 


Substituting equation (7) into equation (10) yields 


2M 

m=l 
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Hence, compatibility is satisfied exactly and the volume integral vanishes 
exactly if satisfies 


= 0 


( 12 ) 


Equation (12) is a well-understood partial differential equation, and con- 
sequently the calculation of q^^^ for practical classes of problems is straight- 
forward. In the next section, values of q^^ are selected for liquids con- 
tained in arbitrarily shaped shell-of-revolution tanks. 

As a consequence of equation (12) , equation (9) becomes 
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6b, 
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dS, 


0 


(13) 


Equation (13) is recognized as an approximate statement of the interface con- 
tinuity between liquid and shell. 

By selecting each series term on the basis of equation (12) , each term 
captures a facet of the liquid behavior observed in classical liquid-shell 
vibration problems, such as periodically supported cylindrical shells contain- 
ing incompressible liquid. (See also the discussion in ref. 7.) In other 
words, each term of the expansion may be viewed as being the exact modal pres- 
sure in a mode of some classical liquid-shell problem. The superposition of 
these terms is forced through equation (13) to approximately satisfy the liquid- 
shell interface condition of tank geometries radically different from the 
classical problem. 


Liquids in Shells of Revolu tion 
General Solution of Laplace * s Equation 

In many practical problems, the liquid is a body of revolution contained 
in a shell of revolution, as shown in figure 1. For such cases, each term of 
the series may be further expanded into trigonometric harmonics about the liquid 
circumference. The deformations of the shell may also be so expanded, and the 
expansion is uncoupled for each trigonometric harmonic. Then, all possible 
solutions of equation (12) , when expressed in cylindrical coordinates, x, r, 
and 9, are the superposition of 
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sh (x^|Xj ) 



( 17 ) 

(X^ > 0) 
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where and are modified and unmodified Bessel functions of 

the first and second kinds, and for the present, it suffices to consider 
arbitrary. For nontoroidal shell-of-re volution tanks, the Bessel functions of 
the second kind, and K^, and r"^ must be dropped since they are singu- 

lar at the axis of revolution (r = 0) . 

Of course, the general solutions of equation (12) need not be expressed in 
cylindrical coordinates, that is, as cylindrical ha 2 rmonics. For example, in 
tanks whose geometry more closely resembles spheres than cylinders, use of 
spherical coordinates should cause the pressure series expansion to converge 
more rapidly. Furthermore, if the liquid is outside of the tank;, it is best to 
use spherical coordinates and hence spherical harmonics for gjnn- However, as 
is shown in the applications sections of this paper, cylindrical harmonics for 
^mn Provide accurate results for spherically shaped tanks with only a few 
terms being used. Further, the cylindrical formulation has the decided advan- 
tage that computer routines for Bessel functions are more widely available than 
routines for the associated Legendre polynomials required if spherical harmonics 
are used. 


Satisfaction of Interface Continuity Relation (Eq. (13) ) 
Substitution of equations (14) to (17) in equation (8) yields 
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(V ^ 0) 


( 20 ) 
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(X^ = 0) 


Again, neglecting the effects due to ullage volume, p — 0 on the free 
surface (x = 0) . Enforcing this condition yields from equation (18) 

M 

m=l 


Inasmuch as equation (21) must hold for all r 




m-HM 


0) 


2 


Substituting equation (22) into equation (7) yields 
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m=l 


g^m 


^ ^ [gmn ^^2 ^ni+M)^m 
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Substituting equations (22) and (23) into equation (6) and using the Tong 
hypothesis of equation (A18) yields the oscillatory pressure. 
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Finally, substituting equation (22) into equation (13) yields 
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(24) 


(25) 
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Since equation (25) holds for all variations of 
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ds 
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( 26 ) 


The following relations nondimensionalize equation (26) : 


r = r/h; x 

= x/h; N = N/h; 

w = w/h; R = R/h; 

s = s/h 

K = 

V = = 

n 

= g/ho)^ 


(27) 


where h is the liquid height, R is the value of r at the shell wall, and 
s is a meridional coordinate as shown in figure 1. Equation (26) now becomes 
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R(s) ds(x) = 0 


(28) 


Following reference 7, w is discretized by assuming a linear interpolation 
of w between liquid nodal stations on the wetted shell wall. Then 

equation (28) may be expressed in matrix notation as 


A - TlA^ - 
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where 1 = k = M and 1 = m = M. The elements of the matrix Q are derived as 
in reference 7 by discretizing the wetted shell wall as 


2ki - (l - -^k,!;,) J__ 
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i+1 _ / X “ Xj^ \ 

R 1 - z ^r^kn <^Vcos \p 
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i-1 \^i - ^i-1, 


X - _ , 

= J^'Jkn dx/cos ip 


where the Kronecker delta is 


(i 7^ j) 


(i = j) 


and i|j is the angle formed by N and the positive r axis as shown in fig- 
ure 1. Also, the notation {w} and {b} represent column vectors of discrete 
normal shell displacements and series amplitude coefficients, respectively. 

y\ y\ 

The elements of and Q 2 are found from equation (34) by replacing 

/- ^‘^kn^ , (- Sgk+M,n\ 

•^kn ^^k^k+M,n - >^x ^ j' respectively. In appendix I 

/N. /\ 

A, A^, A 2 , and A^ are shown to be symmetric matrices. 
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Some insight into the behavior of equation (29) may be revealed by con- 
sidering two limiting cases. In the first case when (jO^ 0 (hydrostatics) , 
equation (29) reduces to 



Substituting equation (36) into equation (24) gives the vector of pressures at 
each liquid nodal station on the wetted shell wall. 



(37) 


where is a gravitational stiffness matrix defined as 


[kJ = gh[B’TA3^Q2] 


(38) 


The elements of the matrix B and the diagonal matrix T are 


- ^%n+M /- - 


B. = 
rm m 3 n 


(xi,Ri) 


(39) 


"ii = 


(i j) 
(i = j) 


(40) 


with Xj^ and being the values of x and r on the shell wall at the ith 

liquid nodal station as discussed in reference 7. The negative sign in equa- 
tion (37) appears because positive pressure indicates compression. 

As a second limiting case, consider the high-frequency oscillations when 
(ij2 co_ Equation (29) then reduces to 


{b> = [a-1q]{w} 


(41) 


Substituting equation (41) into equation (24) gives the vector of pressures at 
each liquid nodal station on the wetted shell wall. 


{p} = 


(42) 
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where M£ is a liquid apparent mass matrix defined as in reference 7 to be 


[%] = [b"a 1q] 


( 43 ) 


where the elements of B** are 


bV = q f X . ,R . ) 

im ^mnV i' 1 / 


(44) 


In general, equation (29) may be solved for the amplitude coefficients as 
follows : 


{b} = [a - nAi - n^A2 - n^A3]“‘ [q - nQi - n^Q2] (w} 


(45) 


provided P = ri* where n* is an eigenvalue of the transcendental equation. 


[A]{b> = + n*^A2 + n*^A3]{b} 


(46) 


Most of the eigenvalues of equation (46) provide the sloshing frequencies of the 
liquid in a rigid tank, that is, for w = 0. Their evaluation is discussed in 
the applications section of this paper. However, equation (46) by itself can 
result in a set of extraneous eigenvalues. This can occur because the varia- 
tional principle used herein and derived in appendix A requires that 6p not 
be identically zero for all x and r. Usually, this poses no problem; how- 
ever, since 6p depends explicitly on the eigenvalue in this formulation, then 
eigenvalues which make 6p vanish identically may exist. In other words, the 
roots of the second bracketed term in equation (28) show up as erroneous eigen- 
values of equation (46) . Such roots occur only when 


n ^ 2 1 

n„ = 1 and T] = - - 

O’ c 
^m^m 


(m = 1 , 2 , . . . ,M) 


(47) 


Hence, erroneous eigenvalues can only appear when gravitational effects are 
included. These eigenvalues may be readily recognized from equation (47) under 
the rare circumstances in which they can occur. 

Substituting equation (45) into equation (24) gives the pressures on the 
wetted shell wall as 


{p} = Pq03^[x]{w} 


(48) 
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where 


[x] = [b][a - - t )^^2 - ^[2 - nQi - 


( 49 ) 


and the elements of w are 


®im gjijj '^mn(^i'%) ’^'^m+M,n(^i'^i)J 


The forces on the shell wall may also be calculated as 


{f} = [t] ^{p} = ^x]{w} 


( 50 ) 


where as in reference 7, 


[tJ = 


a Sj^ 


R 2 As2 


( 51 ) 


Ry ASy 


As a consequence of the Maxwell-Betti reciprocal theorem, the matrix 
product T"^X is symmetric. 


Convergence of Pressure Series 

The convergence properties of the pressure series is a critical subject. 

A general proof showing its convergence in a coupled hydroelastic problem is 
not available and may not be possible. However, it is possible to demonstrate 
its rapid convergence to accurate results through the use of examples. Such 
examples are presented in the applications section of this paper. 
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Hydroelastic Modes 

Hydroelastic inodes are those inodes in which there is significant coupling 
of the liquid and shell motions. In contrast , slosh modes involve little or no 
shell motion- Previous authors have reported that gravity has little effect on 
the hydroelastic modes. (See ref. 1^ ch. 9, and ref. 11^ for example) . In 
other words r| « 1. Consequently/ when searching for the hydroelastic modeS/ 
higher order terms in Tl may be dropped in equation (46) , resulting in 

[x] = [mjQ + n[r] (52) 

where 

[r] = [ba"^] [a^A-Iq - A"1qJ (53) 


Equations (48) and (50) may now be rewritten as 
{p} = 

{f} = + ^[T"^r]){w} (55) 


Equation (54) or (55) provides the wall-deformation-dependent loading terms in 
the shell equations. 


MODIFIED SHELL EQUATIONS 
Discretized Foirm 


If the shell equations have been discretized, as they would be for a 
finite-element or finite-difference shell solution, the liquid loading is 
readily accounted for in the shell equations as. 
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^w 


0 = 


T • 

K K 
UW I ww 





- 0)" 


M, 


uu 


I 

0 M + p T -^Mfj 
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(56) 
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where the K*s refer to stiffness and the M*s to mass and the dimensionless 
vector U includes tangential shell displacements as well as those radial ones 
not on the wetted shell surface; in deriving equation (56) , the approximate 
hydroelastic liquid loading of equation (55) has been used. The eigenvalues 
and eigenvectors can be extracted from equation (56) by the standard 
procedures. 


Differential Equation Form 


In some existing computer programs for shell vibration analysis, the dif- 
ferential equations are solved by numerical integration. This is the case for 
the numerical analysis of references 8 to 10. This analysis and the associated 
computer program have been modified to incorporate the shell-wall-deformation- 
dependent liquid pressure loading developed in the previous section. The 
results presented in the applications section of this paper were generated from 
this modified shell-of-revolution program. 

The program uses a generalized form of Novozhilov*s shell equations which 
include the nonlinear case of moderate rotations. These equations are trans- 
foinned into a set of eight first-order differential equations in eight basic 
force and displacement shell variables and are solved using a Zarghamee version 
of the forward integration method. 

Among other features, the program has branched shell capability and can 
account for effects of prestressing of the shell walls. The branched shell 
capability makes it possible to treat assemblies of tanks often found in 
practice. Prestress effects in shell walls, which effectively stiffen the 
tank, can result from initial pressurization of the tank or the weight of the 
liquid under gravity or other longitudinal acceleration. (Previous work for 
cylindrical shell structures, such as that in ref. 12, has shown that prestress 
effects can be important.) 

The shell differential equations as presented in reference 8 are 


dyn 

ds 


+ a^(s) - b^(s) Zj, = F^(s) 


(57) 


r 


dz^ 


ds 


Yn ^n^^^ ^n ^ 


( 58 ) 


where s is a meridional coordinate. The force vector 
vector Zj^, and coefficient matrices a, b, c, and d 
reference. In addition, the surface load vector Fj^(s) 


y^, displacement 
are defined in the 
is given as 
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and 


X^^> = 0)2103 
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s(l) 


Pn(Sr?) W(C) R(5) 


( 60 ) 


where rUg is the shell mass per unit of shell area and is the nth 

trigonometric circumferential component of the pressure at some meridional loca- 
tion, S/ due to a unit of displaced liquid volume at meridional location, 

To evaluate , w(s) in equation (28)^is set equal to a Dirac delta 

function which vanishes everywhere except at s = f, where it goes to infinity. 
Performing the integrations in equation (28) yield 
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Solving for {b} in equation (61) and substituting into equation (24) result 
in 
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Since A, A^, A 2 , and A 3 are symmetric matrices, it follows from linear 

algebra that Pj^(s,5) is symmetric, that is, Pj^(s,5) = Pj^(C/S). (See, for 
example, ref. 13.) 


APPLICATIONS OF THEORY 

The purpose of this section is to validate the pressure-series theory 
developed in this paper and to demonstrate its rate of convergence through 
appropriate examples. Two sets of examples are considered. The first deals 
with the sloshing of incompressible liquids in rigid cylindrical and hemisphe- 
rical tanks. For the cylindrical tank, it is shown that the present theory 
yields the exact results; while for the hemispherical tank, correlation with 
other approximate theories provides validation. The second set of examples 
deals with the hydroelastic vibrations of incompressible liquids in a flexible 
propellant tank and in a flexible hemisphere. Validation is provided through 
correlation with NASTRAN analysis and experiment. For the hydroelastic results, 
the modified shell-of-revolution program discussed in the section entitled 
"Modified Shell Equations" was used. 


Sloshing Motion of Incompressible Liquids in Rigid Tanks 
General Considerations 


From equation 
given by the roots 

Det Fa - 


(46) , the sloshing frequencies for incompressible liquids are 
of 

* 2 '^ *3'^ n 

n - n A 2 - n A 3 I = o (63) 


/N /N /S 

where the elements of A, A]_, A 2 , and A 3 depend on the values of Aj^^. These 
values are arbitrary; however, judicious choices provide rapid convergence of 
the solution. An examination of some practical cases illustrates the manner in 
which may be selected. 
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Cylindrical Tanks 


For cylindrical tanks, equation (63) greatly simplifies if 
(m = 1,2,.,.,M) are the negative roots of 


d? 


-(?) = 0 


where ^ Then from equations (30) to (33) and (14) to (17) , the 

integrals need only be carried out at x = 1 where n^^ = 1 to give 


pR(l) 

J ('lm+M,n'3kn)^=l^ 


^km ^m 


where 


l^k (1) (1) <kn 


(k m) 
(k = m) 


^l,km \ -0 r 


1 


(k 7 ^ m) 
(k = m) 


0 

^2, km ~ \ _2_ f = > (r-) 

\akCk yr (1) (1) <kn 


(k ^ m) 
(k = m) 


^3 , km 


^k^ktk^^ (1)] <kn 


(k m) 
(k = m) 


K 


= I 


R(l) 


kn-J ^ 

0 


Substituting equations (65) to (68) into equation (63) gives 


(1) - n"a,, y],"' (Dj Li ‘ 


*— ( s) 


.)] [i - n*^5j] = 0 


(64) 


(65) 


(66) 


(67) 


( 68 ) 


(69) 


(70) 
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Hence, setting the first bracketed term of the left side to zero gives 


n* 


coth 


( 71 ) 


which are the exact sloshing eigenvalues for cylindrical tanks. (See ref. 15, 
for example.) The roots of the second bracketed term are, from equation (47), 
associated with the vanishing of 6p and are thus ignored . 


Hemispherical Tanks 

Since the selection of on the basis of equation (64) led to exact 

results for cylindrical tanks, it is reasonable to use this same selection basis 
for geometries of other tanks and allow the superposition of the series terms 
to satisfy the interface condition approximated by equation (13) . To verify the 
validity of this approach, consider a hemispherical rigid tank of radius a and 
liquid height, h = a. Numerical solutions for M = 2, 9, and 10 are given in 
table I along with corresponding solutions from the literature as compiled in 
reference 6. It is seen that the series solution agrees reasonably well with 
published results and converges rapidly. 


Hydroelastic Vibrations of Incompressible Liquids 


in Elastic Tanks 


General Considerations 

Previous investigations have concluded that gravity effects are almost 
always negligible in the hydroelastic vibration modes (ref, 1, ch. 9, and 
ref. 12, for example) - Hence in this section, these effects are neglected 
(i.e., 0=0). From equation (62), the pressure loading on the shell is 

then 



^‘3ln 

T 


Pn(s,C) = Po^'^ ( 

^2n 


'32n 






where the matrix A, as defined in equation (30) , depends on which in 

turn depends on the selected values of (see eqs. (14) to (17)). The 
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terms contain modified or unmodified Bessel functions depending on 

whether is positive or negative, respectively, and contain powers of r 

if Ajj, = 0. 

Following the line of reasoning in reference 7, one selects the positive 
set of as 



and the negative set is selected to satisfy equation (64) . Subsequently, the 
number of positive, negative, and zero values of X^ selected are denoted by 
Ml, M 2 / and M 3 , respectively, whose sum is M. 


Flexible Cylindrical Tanks 

The motivation for using equations (64) and (73) to select A^ is that 
these values lead to a pressure series which provides the exact solution for 
cylindrical tanks. (Compare pressure series with that of ref. 15.) In the next 
section, it is shown that even if the tank geometry is far from cylindrical, the 
superposition of these terms still yields accurate results when applied to the 
interface continuity condition of equation (13) . 


Flexible Propellant Tank and Comparison With Test 

The propellant tank shown in figure 2 has been chosen for the purposes of 
this paper because it comprises a variety of geometries usually encountered in 
aerospace applications- The tank is composed of four distinct sections: a 

nearly spherical upper dome, a conical section, a cylindrical section, and a 
nearly ellipsoidal lower dome. A flange at the juncture of the lower dome and 
the cylindrical section of the tank is clamped to a heavy steel fixture shown 
in figure 2; thus a fully clamped condition was assumed in the analyses. The 
1 / 8 - scale tank approximates a liquid oxygen tank proposed for the space shuttle. 
Its specifications are given in appendix D. 

Nonaxi symmetric and axisymmetric vibrations of the tank were investigated 
both analytically and experimentally. Both the present analysis and the NASTRAN 
hydroelastic analysis were used. The NASTRAN modeling is discussed in appen- 
dix C and the experimental procedure is discussed in appendix D. 

Nonaxi symmetric modes .- Experimental and analytical results for the non- 
axisymmetric modes of the propellant tank are compared in figures 3, 4, and 5. 
Natural frequencies are plotted as a function of the circumferential wave 
number n in figures 3(a), 4(a), and 5(a) for fill conditions of empty, three- 
quarters full, and full. Representative meridional mode shapes for a few select 
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values of n are illustrated in figure 3(b) for 1 and 2 half-waves along the 
shell meridian (m = 1 , 2) and in figures 4(b) and 5(b) for 1, 2 , and 3 half- 
waves along the shell meridian (m = 1, 2 , 3). 

Natural frequencies predicted by the present analysis are in excellent 
agreement with the experiment at all fill conditions and, in general, are some- 
what more accurate than the NASTRAN results- Both the present analysis and the 
NASTRAN analysis included the effects of hydrostatic pressure induced on the 
shell wall from the weight of the liquid- The hydrostatic pressure produces a 
stiffening of the shell which increases with n. This effect is discussed more 
fully in a subsequent section of the paper. 

For the empty tank, the curves of frequency versus the number of circum- 
ferential waves n for m = 1 and m = 2 cross at n = 12, whereas for an 
empty cylindrical tank, these curves, as would be expected, converge at higher 
values of n. The reason for this atypical behavior exhibited in figure 3(a) 
becomes apparent in figure 3(b)- At n = 2, the lowest meridional mode is 
clearly m = 1 with nearly equal response on the cylindrical and conical sec- 
tions of the tank. As n increases, however, the response of the cylindrical 
section decreases until at n = 16, virtually all of the response is in the 
conical section. Similarly, at m = 2, the response of the conical section 
decreases as n increases until at n = 16 all the response is in the cylin- 
drical section- Thus at the higher circumferential wave numbers, the cylindri- 
cal and conical sections of the empty propellant tank behave essentially as 
separate entities. 

The meridional mode shapes from the present analysis, which are virtually 
identical to those generated by NASTRAN, are shown in figures 3(b), 4(b), 
and 5(b)- Qualitative experimental mode shapes were in agreement with the 
analytical results - 

Axi symmetric modes .- Variation of the analytical and experimental fre- 
quencies for the first three axisymmetric (n - 0) vibration modes of the pro- 
pellant tank are presented in figure 6(a) as a function of the liquid depth 
relative to the tank length (h/1) . The agreement of both analyses with experi- 
ment is generally good. The undulating nature of the curves has been ^observed 
by others (ref. 16) and is thought to be a consequence of changes in the 
meridional mode shape . 

The first three analytical meridional mode shapes are presented in fig- 
ure 6(b), for several liquid depths. For clarity, only the mode shapes from 
the present analysis are shown. NASTRAN results are only slightly different 
except at the apex of the lower dome. As noted in appendix D, experimental 
axisymmetric nodal patterns could be detected only on the lower dome. These 
experimental lower dome nodal patterns are identified in figure 6(a) by the 
different test-point symbols and correlate well with the analytical lower dome 
mode shapes of figure 6 (b) . 

Hydroelastic pressure effects in 1/8-scale tank .- As noted previously, the 
analytical results presented in figures 4(a), 5(a), and 6(a) include the effect 
of hydrostatic pressure. The magnitude of this effect on the nonaxi symmetric 
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shell frequencies of the propellant tank model may be observed in figure 7 by 
comparing the results of the present analysis including hydrostatic pressure 
(middle curve) with results from the same analysis without hydrostatic pressure 
(lower curve) . Although the effect of hydrostatic pressure becomes substantial 
for n > 6, its effect is minimal at the lower values of n. If only the lower 
modes are of interest, as is often the case, the mathematical model may be 
simplified by neglecting the hydrostatic pressure with little loss of accuracy. 

Hydrostatic pressure effects in equivalent full-scale tank .- Replica 
scaled models are commonly utilized to study the dynamic behavior of a structure 
in one of two ways- If the model includes sufficient structural detail, dynamic 
analysis may be foregone and full-scale behavior extrapolated from subscale test 
results- More commonly, dynamic tests of a simplified structural model are used 
to validate a mathematical model which in turn is applied to the full-scale 
structure- Caution must be exercised in using either approach to study the 
dynamic behavior of full-scale propellant tanks - 

In general, the natural frequencies of a replica model are inversely pro- 
portional to the scale factor if the internal pressure of both model and full- 
scale counterpart are identical at corresponding locations- This requisite 
causes no difficulties in testing empty tanks but it would be necessary to test 
a 1/8-scale replica model containing liquid in a simulated 8g gravity field if 
the hydrostatic pressure of a full-scale tank at Ig is to be simulated. 

Thus, substantial errors may result if hydrostatic pressure is neglected 
in applying the present analysis to an equivalent full-scale tank containing 
liquid. The magnitude of the effect of hydrostatic pressure for an equivalent 
full-scale tank filled with a liquid is shown by the difference between the 
upper and lower curves in figure 7. (Due to the scaling factor, all the fre- 
quencies indicated in this figure are 8 times those of a full-scale tank.) 

The upper curve results from the present analysis of the subject propellant 
tank (1/8-scale) in an 8g gravity field. As may be seen, the effect of hydro- 
static pressure is much greater for the full-scale tank than for the 1/8-scale 
model. For n = 6, neglecting hydrostatic pressure effects results in a fre- 
quency prediction error of only 15 percent for the 1/8-scale model subjected to 
Ig but nearly 100 percent error for an equivalent full-scale tank at Ig . These 
scale effects of hydrostatic pressure would become even more pronounced if the 
longitudinal acceleration during flight was taken into account . 

Free surface motion .- Figure 8 illustrates the free surface mode shapes 
for different numbers of circumferential waves n for m = 1, 2, and 3 merid- 
ional half-waves in the full and three-quarters full propellant tank. As indi- 
cated in figure 8, the entire free surface is moving when n = 0; however as 
n increases, the motion becomes localized to the shell wall with little motion 
occurring near the center line of the tank. This phenomenon is known to occur 
in practice and is readily accounted for in the present analysis by the modified 
Bessel functions in the series expansion. This helps to explain why the series 
expansion is so accurate for this type of problem. 

It is also interesting to observe that the surface motion undergoes little 
change with meridional wave number. Consequently, mode shapes must be charac- 
terized by the motion of the shell wall. 
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Convergence of pressure series solution .- For the full propellant tank/ 
analyses were performed using several combinations of M 2 / and M 3 in 

order to study the significance of each type of series term in such problems. 

The analyses were performed for the three lowest frequencies at n = 6 , 
m = 1 , 2 , and 3 ; frequency results are presented as ratios to the corresponding 
converged values in table II. These results show that with M-|_ = 5, converged 
frequencies are obtained irrespective of the values of M 2 and M 3 . (As 
before, M]_, M 2 , and M 3 are the number of user-selected positive, negative, 

and zero where the positive A^^ are given by equation (73) and the nega- 

tive A^i by equation (64).) Hence, only the positive set of Aj^ need be 
selected, although addition of unnecessary terms did not deteriorate the 
results. Though not shown in table II, this conclusion has generally been true 
for all n > 0 (nonaxi symmetric modes); however, for n = 0 (axisymmetric 
modes), Aj^^ = 0 should be retained (i.e., M 3 = 1 ) . 

In figure 9, the series convergence at n = 6 is further examined for 
increasing values of M^_ with M 2 = M 3 = 0. It is clear from the figure that 
the convergence of the series is very rapid. Consequently, it is concluded that 
although the Ajj^ are selected on the basis of hydroelastic vibrations in cylin- 
drical tanks, the same selection may be safely extended to other geometries . 

Note that for nonaxisymmetric modes, the lower dome has negligible motion; 
hence the validity of the present analysis for such a dome undergoing nonaxi- 
symmetric motion is not confirmed by the propellant tank results- To confirm 
the validity of the present analysis for such a dome, a partially filled elastic 
hemisphere is considered. 


Flexible Hemisphere 

The first two free edge modes of a partially water-filled aluminum hemi- 
sphere of radius a having two circumferential waves (n = 2 ) are shown in fig- 
ure 10 for the full spectrum of fill conditions, 0 ^ h/a ^ 1. Frequencies and 
mode shapes predicted using the present analysis are compared with NASTRAN 
results. For this comparison, the pressure loading on the shell wall, generated 
by NASTRAN, was applied to the structural program (refs. 8 to 10) used by the 
present analysis. Thus, any discrepancy between the two results can be attri- 
buted only to differences between the liquid finite-element analysis employed 
in NASTRAN and the pressure series analysis proposed herein. Since results for 
the propellant tank, previously discussed, indicated that only positive values 
of selected on the basis of equation (73) need be considered, only series 

terms corresponding to these values were used for the hemisphere. The conver- 
gence of the series expansion is indicated in figure 11. As in the case of the 
propellant tank, few terms need be retained. 

In figure 10 the converged solutions of the present analysis for frequency 
and mode shape generally correlate very well with converged NASTRAN-generated 
solutions. This example further indicates that the use of a series containing 
products of modified Bessel functions and trigonometric functions ideally 
applicable to liquids in circular cylindrical tanks can be superimposed to 
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solve other geometries with adequate accuracy. Discrepancies between the 
two solutions in the second mode at low liquid height probably are due to the 
scarcity of liquid finite elements in the NASTFAN model near the bottom of the 
tank. 


CONCLUDING REMARKS 

An analysis is presented of hydroelastic vibrations of elastic tanks par- 
tially filled with incompressible nonviscous liquid in the presence of gravity. 

A complementary energy principle suitable for handling the liquid is derived, 
and the liquid is analyzed by using a series representation of the oscillatory 
pressure. By selecting each term of the series to satisfy the liquid governing 
equation (Laplace's equation), the energy principle reduces to an integral 
equation on the wetted portion of the tank wall. This integral equation 
represents an approximate statement of the liquid/elastic- tank compatibility 
at the wetted tank wall and may be solved numerically, without discretization 
of the liquid volume, to yield the pressure loading on the wall as a function 
of wall motion. 

This approach permits the wall-motion-dependent liquid pressure loading to 
be readily incorporated into existing shell analysis computer programs. This 
procedure was followed to produce a modified shell-of-revolution program for 
hydroelastic vibrations. In the modified program presented herein, the pressure 
series terms (which must satisfy the governing liquid equation) were selected on 
the basis of a circular cylindrical tank geometry; that is, each term is a pro- 
duct of a modified (or unmodified) Bessel function and a trigonometric (or 
hyperbolic) function. The energy principle then allows solution of problems 
with geometries other than cylindrical by superposition of these terms. Appli- 
cation of the modified program to the gravity sloshing or hydroelastic vibra- 
tions of incompressible liquids in tanks of various geometries provides the 
following conclusions: 

1. Since cylindrical functions are used as series terms, exact solutions 
are obtained for cylindrical tanks. 

2. Sloshing frequencies for rigid hemispherical tanks correlated very well 
with published frequencies; and hydroelastic vibrations in an elastic propellant 
tank and in an elastic hemisphere also agreed well with test and finite-element 
(NASTRAN) results- Thus, the procedure can be applied to a wide variety of tank 
geometries . 

3. Accurate results require retaining only a few terms in the series 
expansion. 
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4- Selection of the series terms is straightforward and is readily incor- 
porated as an automatic feature of the modified program. 

In addition, it is found that hydrostatic pressure from the liquid weight can 
cause significant stiffening of the tank wall and must therefore be accounted 
for, especially in full-scale fuel tanks on large space vehicles. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 30, 1979 
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APPENDIX A 


LIQUID COMPLEMENTARY ENERGY PRINCIPLE 

The purpose of this appendix is to formally state and confirm the validity 
of the liquid complementary variational principle used in the analysis presented 
in this paper. This principle is analogous to such principles in structures. 

The liquid pressures (or in structures, the stresses) are varied within a class 
of admissible pressures (of stresses) ; the admissible variations must satisfy 
Euler's equations (or the equations of motion) for perturbed motions about an 
initial state. When the variations are carried out, the complementary principle 
yields a compatibility equation (or equations) and boundary conditions on veloc- 
ities (or displacements) . 

This appendix proceeds by first stating the principle, then establishing 
the admissible class of variations, deriving the appropriate variational comple- 
mentary energies, and finally proving the principle. 


Statement of Principle 

Consider a nonviscous incompressible liquid which . undergoes small irrota- 
tional harmonic oscillations, with the complex form about an initial 

pressurized state and in the presence of a gravitational field. The deformation 
state satisfies the compatibility relation, 


in 


= 0 

and 

3 ( 1 ) 

= w 

3n 


(Al) 


(A2) 


on 


where w is applied. This deformation state makes the variation of the 


total complementary energy, vanish, that is. 


= 0 


(A3) 


provided that the variations satisfy Euler's equation and 


6p = 0 


on S^. The next section of this appendix shows that constraining the varia- 
tions to satisfy Euler's equation is equivalent to their satisfying a linearized 
form of Bernoulli's equation for small perturbations, namely. 


5p = -Pq6(J) + Pog<Su^ 


(A4) 
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Furthermore, 


(A5) 

c c c,g c ' ' 

where T^, and W_ - the complementary kinetic, gravitational, and 

C- ^ f y ^ 

external work energies of the liquid, respectively - are derived in a subse- 
quent section of this appendix. 


Variation of Euler's Equation for Small Perturbed Motions 

Euler's equation which expresses dynamic equilibrium for a nonviscous 
incompressible liquid is (see, for example, ref. 17, p. 95) 




+ p, 


(u • div)u + 


grad P - 


Pq5 


grad x = 0 


(A6) 


where the local Eulerian variables are as follows: p^, the uniform mass density 

of the liquid; U, the velocity; and P the pressure. Considering a small 
perturbation about an initial pressurized state, one can write 




U = u 


P = 


Pq P 


x = 


^o + ^x 




(A7) 


In equations (A7) , subscript o refers to the initial state and it is assumed 
that the liquid is initially at rest. Substituting equations (A7) into equa- 
tion (A6) and neglecting nonlinear terms because the perturbations are assumed 
small yield 


p^u + grad p - p^g grad u^ + (grad - p^g grad x^) = 0 (A8) 


The term in parentheses vanishes since it represents equilibrium in the initial 
state and equation (A8) becomes 

Po^ + grad p - p^g grad = 0 (A9) 
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By introducing the velocity gradient for irrotational flow, as given by equa- 
tion (5) of the main text, into equation (A9) , 

Pq grad ^ + grad p - p^g grad u^^ = 0 


This equation may be integrated spatially to produce a linearized form of 
Bernoulli's equation for perturbed liquid motion, namely. 


Po<i> + P - Po9^x = 0 


(AlO) 


where the constant of integration may be safely incorporated into (p as dis- 
cussed in reference 1, page 15 and elsewhere. The variation of equation (AlO) 
may now be taken to yield equation (A4) . 


Variational Complementary Forms of Kinetic, Gravitational, 

and External Energies 

Complementary kinetic energy .- The well-known complementary kinetic energy 
is given as 


T 


c 



u dv^ 


(All) 


where the integration may be safely carried out over the undeformed volume 
since the liquid is assumed incompressible. 

Substituting equation (5) into (All) yields 


T 


c 



Pq grad (j) • grad (p dV^ 


(A12) 


Taking the variation of both sides of equation (A12) results in 


6t^ = Pq grad (p “ grad 6(p dV^ 


(A13) 
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Green’s theorem then provides 


<STc 



Po 64) dVQ 



pQ grad (j) • N 5 (j) dS^ 


(A14) 


where N represents an outward unit vector normal to the undeformed liquid 
surface S^, 

Complementary gravitational energy ,- The variation of the complementary 
potential energy is found by considering the perturbed state of the liquid and 
varying the gravitational body forces which are given by 


F = P^g grad x 

Taking the variation of F and using equation (A7) gives 
5F = p^g grad (6uj^) 

Hence the variation of complementary gravitational energy is 

6Uc q = \ 6 f • u dVo = \ Pq 9 grad 6ux • u dVQ (A15) 

'^Vo 


Substituting equation (5) into equation (Al5) yields 


6u 


c,g 


j03 


1 


pQg grad 6u^ • grad (p dV^ 


(A16) 


Application of Green's theorem to equation (A16) gives 


6u 


c,g 


^ r 


V^<t> 6ujj 


dV, 



8n 




dSo 


(A17) 
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Some previous authors (e.g., refs. 18 and 19) have used the Tong hypothesis, 
which states that 




n^(u 


303 oN 


(A18) 


where is the cosine of the angle between the local surface normal N and 

the x-axis. Under this hypothesis, the surface integral of equation (A17) 
becomes 



9n 




dS, 



9n 


6 



dS 


o 


Substituting equation (A19) into equation (Al7) gives 


6u 


c,g 



V2fj) 6u^ dVg 



9n 



(A19) 


(A20) 


Thus equation (A20) establishes the gravitational energy variation for a per- 
turbing oscillation about an initial state of an incompressible liquid. 

Complementary external work ,- The variation of the complementary external 
work i s 


6Wc 



w 6p dS^ 


{A2D 


where w is the perturbed displacement at the liquid surface in the N 
direction. Substituting for 6p from equation (A4) , one has 


SWc = 

Use of the Tong 


6Wc = 


-pQ w(-jw 64) + g 6uj^)dSQ 


hypothesis (eq. (A18) ) yields 


jf^Po 


1 


o L 


gn 

6(J) + — S 


0)^ 


m 


dSr 


(A22) 


(A23) 
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Proof of the Variational Principle (Eq. (A3) ) 

Siibstituting equations (A14) , (A20) , and (A23) into equation (A3) results 


in 


0 = 611^ = ^ [po v2(()(-jto 6(}) + g dV^ 


+ Pc 




dS, 


(A24) 


Use of equation (A4) in the volume integrand of equation (A24) gives 


0 = 6n, 


= r v2(f, 


6p dVp + pQ 


r 

Jq \ 9N 


j(JL)W 


S(t> + 


gn 




OJ 


2 


dS. 


(A25) 


Since equation (A25) must be valid for all variations. 


V^cj) = 0 


(A26) 


in and 


3cj) 

= W 

9n 


on S and 
w 


&<i> + 



0 


(A27) 


on which from equation (A4) is satisfied by the constraint 6p = 0 on Sf. 

Hence, the principle is verified. 


34 



APPENDIX B 


PROOF OF SYMMETRY OF MATRICES A, A^, A 2 # AND A 3 

To prove the symmetry of A as defined in equation (32) , it suffices to 
show that 

X 

Inasmuch as for 1 = m = M vanishes on the free surface (x - 0) , the 

integration of equation (Bl) may be carried out over the entire surface, Sq. 
Green's second formula then gives 

r (%n " ^itin = ° 

^o 

But equation (B2) is valid as a consequence of equation (11) ; hence the proof is 
complete . 

/s 

similar proofs can be worked out for A^, A 2 / and A 3 . 


^qmn 

9n 


3qkn 


(1 = m = M; 


1 = k = M) 


(Bl) 
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NASTRAN MODELS OF PROPULSION TANK AND LIQUID 

Finite-element frequencies and mode shapes used for comparison with results 
of the analysis presented in this paper were obtained using mathematical models 
described in this appendix- These models were developed within the NASTRAN 
hydroelastic capability (ref. 4) and modified according to references 2 and 3 to 
permit a symmetric fluid matrix formulation and thereby simplify computational 
procedures. The analysis is limited to incompressible fluids, and computational 
accuracy and efficiency are significantly increased by introducing a trigono- 
metric series representation of circumferentially distributed displacements. 

For the present purposes, the modified NASTRAN analysis was used only for 
the case of a free liquid surface, with gravity effects neglected. The analysis 
is valid also for the more general situation of a pressurized liquid level with 
gravity effects included. 


Essentials of Modified NASTRAN Analysis for This Study 

The theoretical development in reference 2 is based on the complementary 
energy principle applied to systems of discrete variables. The complementary 
kinetic energy of the liquid is the quadratic function 

■'o " I 2 S '•ijSiPj " 

i j 

written in terms of generalized impulsive pressures which are based on the 
pressure impulse definition 


P 



dt 


(C2a) 


or 


P 



(C2b) 


Since the liquid is assumed incompressible, there is no strain energy. The 
quantity Lj_ j , an element of the symmetric inverse mass matrix (termed inertance 
matrix in ref. 2) , is defined as 




9p. 9p. 

1 J 


(C3) 
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The complementary virtual work for the discretized system is expressed by 


6w^ 


22 


Ai^u^ 6pi 


(C4) 


where is a discrete displacement of the liquid bounding surface and 

is the corresponding elemental generalized area as defined in reference 2. 

Equations (Cl) to (C4) are basic to the derivation of the Euler-Lagrange 
equations of motion of the liquid in reference 2. For the incompressible liquid 
these equations are shown therein to reduce to the partitioned matrix equation. 
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where subscripts f, s, and i denote free liquid surface, liquid-structural 
interface (or wetted surface), and liquid interior, respectively. As a conse- 
quence of incompressibility, p^ may be related to fPf.Ps'l'^ so that 



(C6) 


in terms of the reduced inertance matrix 


— 

— 


~ 

— 


— — 1 

^ff 

^fs 


Lff 

^fs 


^fi 

Isf 



_^Sf 

^ss_ 






(C7) 


For liquid displacements (or accelerations) normal to the free liquid surface, 

the off-diagonal area submatrices A*^ and A*^ vanish. 

fs sf 

Dynamic coupling of the liquid with the tank structure is accomplished by 
introducing pressures determined from equation (C6) into a standard set of 
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Structural dynamic equations, which may be written in the corresponding par- 
titioned matrix form 


0 

0 




(C8) 


Since the entire surface is constrained, that is, 

(Pf) = 0 (C9) 

equation (C6) is solved to give 

and 

(Sf) - W tL](“s) 

Substitution of equation (CIO) into equation (C8) results in 

[^s + Ml](u3) + [Kss](^s> = 0 (C12) 

which is the equation of free vibration of a liquid-loaded tank, with the liquid 
loading given by the liquid mass matrix 


For simple harmonic motion, equation (C12) reduces to the familiar eigenvalue 
equation 





(C14) 


which gives the frequencies and mode shapes of a liquid-loaded tank in free 
vibration. 
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Computer Program Implementation 

The NASTRAN hydroelastic computer program of reference 4 is utilized and 
modified to implement the analysis of reference 2. As detailed in reference 3, 
the unsyiranetric matrix formulation in the complex eigenvalue module (Rigid 
Format 7) is altered to compute the symmetric liquid mass matrix of equa- 
tion (C13) . The combination of this matrix with the tank structural mass and 
stiffness matrices, as indicated in equation (C14) , is programmed in the normal 
vibration mode module (Rigid Format 3) . The program modifications were made in 
accordance with the NASTRAN Direct Matrix Abstration Program (DMAP) for altering 
rigid formats. 


Harmonic repre sentations . - The liquid geometry is assumed to be axisym- 
metric, so that pressure distributions can be represented by a trigonometric 
series, such as (e.g., in ref. 2) 




Po(ri,Zi) 


n 

cos k0j_ + (CIS) 


k=l 


in terms of circumferential harmonics (n = 0,1,2,...), where pj^ and are 

generalized pressure amplitudes of symmetric and antisymmetric distributions, 
respectively. The condition of uniform pressure is given only for n = 0 in 
equation (CIS) . Harmonic pressure distributions such as equation (CIS) are pro- 
grammed in reference 4 and unchanged in reference 2. 


A significant modification introduced in reference 2 is the trigonometric 
series representation of structural grid displacements consistent with equa- 
tion (CIS) . Relationships between discrete-displacement and harmonic- 
displacement amplitudes are listed therein' for cylindrical and spherical coordi- 
nate systems. All rigid-body motions are included in these relations. The 
transformations from discrete to harmonic displacements (in ref. 2) are made by 
including multipoint constraints (MPC) of the form 


('^g) = [°gh]K) <ci6) 

in the input data, where discrete grid displacement vector, 

is the harmonic-displacement vector, and transformation matrix. 

Elements of trigonometric functions evaluated at discrete-variable 

grid locations of the tank structure. Harmonic displacements are oriented to a 
set of fictitious grid points distributed along the tank meridian and distinct 
from physical grid points. At each meridian in the physical grid system, physi- 
cal grid points around the circumference are related to the so-called harmonic 
grid points by equation (C16) . 

This transformation not only makes displacement and pressure representa- 
tions consistent but can also result in large reductions in the number of 
degrees of freedom for hydroelastic dynamic analysis. For all six degrees of 
freedom at each physical grid point, the maximum number of degrees of freedom in 
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the physical system is 6 x J x K, where J denotes the number of meridians and 
K the number of circumferential grid points at each meridian. This triple 
product is the dimension of (i.e., momber of rows) in equation (C16) . 

If N denotes the number of harmonics (or circumferential wave numbers 
n = Orl/2,...) included in the analysis, the dimension of {^^h^ is 6 x n, and 
since N is usually much smaller than either J or K, equation (C16) can 
result in a radical reduction in problem size, namely N/JK. Still further 
reduction is possible by omitting generalized rotation and tangential degrees 
of freedom in a Guyan reduction, with harmonic grid points referenced in 
NASTRAN OMIT instructions. 

Computations of liquid pressures and masses .- Liquid finite elements in 
NASTRAN are rectangular and triangular solids of revolution interconnected by 
concentric rings comprising the liquid grid system. The liquid elements 
(CFLUIDi, i=2,3,4) and their ring boundaries (RINGFL specifications) are analo- 
gous to structural finite elements and grid points, respectively. Liquid incom- 
pressibility is represented by removal of internal pressure degrees of freedom 

in eq. (C5)j by means of a special OMIT instruction in which the con- 
centric rings are referenced by a 7-digit integer numbering system. The surface 
pressure constraints in equation (C9) are represented by special single-point 
constraints (SPC) , also referencing the 7-digit ring identification. The wetted 
tank surface is represented by GRIDB instructions which relate structural grid 
points (GRID) there with liquid-element boundaries. Since liquid pressures 
calculated by equation (CIO) depend primarily on outward normal (or radial) 
accelerations, rotations and tangential displacements are constrained (by SPC 
instructions) , and MPC relations between physical and harmonic degrees of free- 
dom (eq. (C16) ) are reduced by a factor of 6. 

The symmetric fluid mass matrix given by equation (C13) is calculated by 
modifying the program steps of Rigid Format 7 so that the complex eigenvalue 
routine (CEAD) is not used. These modifications are specified in reference 3 by 
DMAP ALTER instructions for matrix printout, partitioning, and matrix inversion 
to form the liquid inertance, surface area, free surface accelerations (or dis- 
placements) , pressure, and mass matrices of equations (CIO), (Cll) , and (C13) . 
Along with these types of instructions are direct matrix input (DMI) specifica- 
tions in the input data in the form of partitioning vectors to identify pressure 
and displacement degrees of freedom. The OUTPUT! DMAP instruction causes the 
liquid pressure and mass matrices to be punched on cards in DMI format. 

Vibration modes of liquid-loaded tank .- Liquid pressures and masses, plus 
free surface motions (eq. (Cll)), as desired, output in Rigid Format 7 are 
input as data in Rigid Format 3 to obtain normal free vibration modes. Liquid- 
structural coupling is directed by DMAP ALTER instructions as shown in refer- 
ence 3 prior to execution of the real eigenvalue routine (READ) . The tank 
structure is modeled by flat plate elements (CQUAD2) interconnected by grid 
points. Harmonic grid points are included, as in the liquid model, and the MPC 
relations include harmonic representations of all six degrees of freedom (three 
displacements and three rotations) at each physical grid point. Rotations and 
tangential displacements, constrained by SPC instructions in Rigid Format 7, are 
removed by OMIT instructions in Rigid Format 3. 
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NASTRAN Models 

NASTRAN finite-element models of the liquid-loaded elastic propellant tank 
of figure 2 are illustrated in figures 12 and 13. The model in figure 12 is 
similar to that used in references 2 and 3 except for the liquid element bound- 
aries in figure 12(a). A finer grid model is shown in figure 13. Free surface 
motions^ given by equation (Cll) , were not sought. Physical and hartnonic grid 
points in the upper and lower domes and in the conical section are oriented in 
spherical coordinates, and grids in the cylindrical section are in cylindrical 
coordinates. The closer spacing between concentric rings near the tank wall is 
intended to simulate more accurately large pressure gradients there, particu- 
larly for higher circumferential wave numbers (or harTuonics) - The meridional 
liquid grid locations match the meridional structural grid locations along the 
wetted surface. 

Liquid models .- The coarse liquid grid system in figure 12(a) consists of 
15 meridional stations with 4 liquid-element boundaries at each meridian except 
3 near the tank bottom which have fewer boundaries. The model has a total of 
1386 possible pressure and displacement degrees of freedom. The free surface 
pressures are constrained by 16 liquid SPC statements, and all but the radial 
displacements are constrained by 987 structural SPC*s. Incompressibility 
accounts for another 132 pressure degrees of freedom being removed by liquid 
OMIT statements. Harmonic representations of radial displacements through MPC 
relations account for another 135 degrees of freedom. 

There are 116 degrees of freedom remaining, of which 56 are pressures and 
60 are displacements. Four harmonics (N = 4) are included in this mode, 
n = 0 to 3, so that (J-l)N = 56 and JN = 60, where J in this case is 15. 
These integer products determine the matrix sizes and partitions for liquid 
pressures and masses on the DMI specifications for the vector P9 in the DMAP 
ALTERS of reference 3. The pressure matrix (PDU2) in reference 3 computed for 
this model is rectangular with 56 rows and 60 columns, and the symmetric liquid 
mass matrix (MFLD) is 60 x 60. 

The fine liquid grid system in figure 13(a) consists of 28 meridional sta- 
tions with 6 liquid-element boundaries at each meridian except 2 meridians 
closest to the tank bottom, where 4 and 5 rings bound 5 liquid elements. The 
maximum number of degrees of freedom for this model is 3537. There are 2818 SPC 
statements, of which 6 are fluid SPC * s constraining free surface pressures. 
Liquid OMIT statements account for 132 degrees of freedom and MPC relations for 
another 532. In the remaining 55 degrees of freedom, there are 27 pressures and 
28 displacements. With one harmonic at a time (N = 1) in this model, the number 
of displacements is the same as the number of meridians (J = 28, N = 1) , and 
the number of pressures is one less. Thus, the liquid pressure matrix is 
27 X 28 and the liquid mass matrix 28 x 28. 

For partially filled tanks, the number of pressures is reduced, but not the 
number of displacements. The reduced number of pressures is obtained by sub- 
tracting more than one meridian from J for the full tank; that is, J-1 
becomes J-2, J-3, J-5, and so forth for various reduced levels of the liquid. 

Pressures above the liquid surface are constrained by liquid SPC statements. 
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Structural models .- The coarse structural model of figure 12(b) has 9 phys- 
ical grid points around half the circumference at each meridian/ with a 22.5^ 

interval between circumferential grid points. The top and bottom domes are 

modeled by triangular plate elements (CTRIA2) which connect the apex grid points 
to the grid points at the full tank level and at the lowest meridian. The maxi- 
mum number of degrees of freedom is 1182, including 60 harmonic displacements. 
Constraints (SPC) are applied to 72 degrees of freedom, of which 48 satisfy 
boundary conditions for motions associated with n = 0 that are symmetric with 
respect to the pitch plane of symmetry. The remaining 24 constrained degrees of 
freedom simulate the tank mounting in the test program (appendix D) , that is, 
the clamped ring modeled by beam elements (GEAR) , at the base of the cylindrical 
section. Application of the Guyan reduction accounts for 240 omitted degrees of 
freedom (on OMIT statements) , of which 28 are associated with harmonic n = 0 

and the rest with harmonics n = 1, 2, and 3. There are 810 MPC relations to 

implement the harmonic representation of all 6 physical degrees of freedom at 
each grid point. The remaining degrees of freedom number 60. Although har- 
monics n = 1, 2, and 3 are included in the calculations, vibration modes of 
tank and fluid are reported only for variations in tank fullness associated with 
the axisymmetric mode n = 0, for which the coarse liquid and structural grids 
illustrated in figure 12 are considered adequate. 

The fine structural model of figure 13(b) has 19 physical grid points 
around a quarter of the circumference at each of 30 meridians, with a 5^ inter- 
val between circumferential grid points. Two of these meridians lie above the 
full tank level, and apex grid points are connected to the highest and lowest 
meridians by triangular plate elements. There are a total of 3612 physical and 
harmonic degrees of freedom which reduce to 28 harmonic degrees of freedom. 
Constraints specifying symmetric boundary conditions at the apexes plus clamped 
ring constraints total 12. The Guyan reduction is applied to the apexes, con- 
sistent with symmetric boundary conditions, for all harmonic degrees of freedom 
at the two meridians above the full tank level and for all but the harmonic 
radial degree of freedom over the wetted tank surface. The number of omitted 
degrees of freedom totals 152. Harmonic representation accounts for 3420 degrees 
of freedom. The fine model was used to calculate liquid-loaded vibration modes 
for nonaxisymmetric modes ranging from n = 2 to 16 for both three-quarters 
full and full tanks. 

Hydrostatic pressure loading .- The effect of hydrostatic pressure was 
introduced into the NASTRAN finite-element model by determining an additional 
stiffness matrix equivalent to the matrix product in equa- 

tion (56) . Static pressures at a given fluid depth were calculated by the pro- 
duct of the depth by the unit weight of liquid and were uniformly distributed 
over quadrilateral plate elements forming the tank boundary at that depth. 

These pressure loads were computed for every row of plate elements from the 
liquid level to the tank bottom and were input (on PL0AD2 statements) to the 
structural model for n = 0. This model was then run in the NASTRAN differ- 
ential stiffness module (Rigid Format 4) which was instructed by DMAP ALTER to 
output CakJ onto magnetic tape. This additional stiffness matrix was then 
read (also by DMAP ALTER) into Rigid Format 3 and added to the structural 
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stiffness matrix |_^ssj • eigenvalue equation (C14) was then solved in the 

usual manner to obtain the natural frequencies and mode shapes of the liquid- 
loaded tank under hydrostatic pressure. The effect of this hydrostatic pres- 
sure loading was determined for full and three-quarters full tanks for nonaxi- 
symme trie modes using the fine-grid model of figure 13. 
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PROPELLANT TANK TESTS 

This appendix contains a brief description of the propellant tank referred 
to in the main text and the equipment and procedure used to obtain the data pre- 
sented therein. 


Propellant Tank 

The propellant tank depicted in figure 2 is a 1/8-scale model of an early 
version of the space shuttle liquid-oxygen tank. The tank is made up of four 
aluminum sections: an upper dome of nearly spherical shape, a conical section, 

a cylindrical section, and a lower dome of nearly ellipsoidal shape. The over- 
all diameter and height of the tank are approximately 1 and 2 meters, respec- 
tively, The wall thicknesses of the various sections are indicated in figure 2. 
A flange at the juncture of the lower dome and cylindrical section of the tank 
is clamped to a heavy steel test fixture shown in figure 2. 


Nonaxi symmetric Tests (n / 0) 

Sinusoidal excitation of the tank was provided by a small servocontrolled 
electrodynamic exciter driving normal to the surface of the cylindrical or 
conical section of the tank and, where applicable, below the liquid surface. 
Maximum force capability of the exciter was 4.4 newtons. A force gage used to 
monitor the input force was installed between the exciter and tank shell. A 
servocontrol and oscillator were employed to maintain a constant sinusoidal 
exciter force. 

The response of the shell was detected by a motorized noncontacting dis- 
placement transducer which traversed a circular track (fig. 2) . The track could 
be positioned at any desired station along the length of the tank. A coincident- 
quadrature analyzer was used to measure the quadrature component of the dis- 
placement (90° out of phase with the input force) . Resonance was determined by 
manually adjusting the exciter frequency at a given level of input force until a 
maximum quadrature component was obtained. The circumferential mode shape was 
then recorded by inputting the quadrature component of the displacement into an 
x-y plotter as the transducer traversed the circumference of the tank. Quali- 
tative longitudinal variation of the mode shapes was obtained by observing the 
Lissajous pattern on an oscilloscope screen as a hand-held velocity transducer 
was moved to several axial stations along the tank. Vibration modes were mea- 
sured for liquid (water) depths of 0, 1.27, and 1.90 meters corresponding nomi- 
nally to empty, three-quarters full, and full conditions. 


Axisymmetric Tests (n = 0) 

For the axisymmetric vibration modes, analyses indicate (fig. 6(b)) that 
the largest vibration amplitudes occur below the liquid level and that at the 
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greater liquid depths, the vibration amplitude of the conical and cylindrical 
sections are of magnitude comparable with that on the lower dome. In the labo- 
ratory, however, the response of the cylindrical and conical tank sections in 
the high-energy axisymmetric modes is completely masked by the large response in 
the neighboring low-energy nonaxi symmetric modes. The problem is further aggra- 
vated for this particular tank by the high modal density of the nonaxisymmetric 
shell modes in the same frequency range as the axisymmetric modes. Axisymmetric 
nodal patterns were identifiable, however, on the lower dome, and in figure 6(a) 
the resonant frequencies of the first three lower-dome n = 0 modes are plotted 
as a function of the liquid depth. Except at the shallower liquid depths, an 
axisymmetric dome response was always coupled to a much larger nonaxisymmetric 
response of the cylindrical and conical sections; thus, most of the data points 
of figure 6(a) are in reality the resonant frequencies of n / 0 shell modes 
in proximity to the resonant frequencies of the n = 0 dome modes. In some 
cases, two or more n = 0 dome responses having the same nodal pattern were 
present in a small frequency band. If the response at one frequency was. appre- 
ciably larger than the others, only one data point was recorded. If the 
response of two closely spaced resonances were of comparable magnitude, both 
frequencies were recorded. 

Many of the axisymmetric modes did not exhibit classical characteristics in 
that where one segment of the lower dome contained only one circumferential node, 
another segment might contain two. Such cases are designated in figure 6(a) by 
the superposition of the appropriate nodal pattern symbols. 

The noncontacting probe and track assembly used for the nonaxisymmetric 
shell modes could not be readily adapted to measure responses on the lower dome. 
The modal patterns on the lower dome were detected either by touch or through 
the use of a movable vacuum-attached accelerometer and Lissajous patterns. 
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TABLE I.- SLOSHING FREQUENCIES IN A RIGID HEMISPHERICAL TANK 




Frequency parameter. 

i/n = to^ 

h/g 


n 

Results 

compiled in 

ref, 6 

Present solution 
(M^ = M 3 = 0 ) 


Finite 

elements 

Method of 
ref . 6 

Integral 

solution 

M 2 = 2 

M 2 = 9 

M 2 = 10 

0 

3.96 

3.92 

3.65 

3.862 

3.791 

3.790 

1 

1.57 

1.68 

1.56 

1.620 

1.594 

1.592 

2 

— 

3.03 

— 

2.925 

2.889 

2.887 

3 

— 

4.34 

— 

4.295 

4.349 

4.359 

4 

— 

5.62 

— 

5.370 

5.411 

5.414 

5 

— 

6.88 

— 

6.445 

6.473 

6.475 


TABLE II,- EFFECTS OF THE THREE TYPES OF SERIES TERMS 

ON THE n = 6 PROPELLANT TANK FREQUENCIES 

Subscript c refers to converged results 
and subscripts 1 , 2 , and 3 refer to 

m = 1 , 2, and 3. 


Ml 

1 

M2 

M3 

•H 

s 

II 

II 

Wl/Wic 

^2/^20 

W3/W3C 

5 

0 

0 

5 

1.00 

1.01 

1.09 

0 

1 

0 

1 

7.04 

7.37 

6.47 

0 

5 

0 

5 

7.06 

6.39 

5.93 

0 

0 

1 

1 

1.74 

5.12 

6.45 

5 

5 

1 

11 

1.00 

1.01 

1 

1.09 
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(b) Analytical meridional mode shapes. 


Figure 3.- Concluded. 
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Figure 4,- Concluded. 
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Figure 6.- Axisymmetric (n = 0) 


Experimental 
lower -dome 
nodal pattern 



Experiment O □ 

Present analysis 

NAS TRAN 





ive liquid height, h/l 
frequency with liquid height. 


vibration modes of propellant tank. 






PERCENT 

ERROR 


m = 1 2 3 



NUMBER OF TERMS IN SERIES. M^ 

Figure 9.- Accuracy and convergence of predicted 
frequencies of first three meridional inodes 
(m = 1,2,3) for full propellant tank. n = 6; 
M2 = M3 = 0. 
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Figure 11.- Accuracy and convergence of predicted 
frequencies of first three meridional modes 
(m = 1,2,3) of full hemisphere. n = 2; 
l2 = M3 = 0. 
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